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Abstract 

We explore some basic entanglement features of multiqubit systems that are relevant for the 
development of algorithms for searching highly entangled states. In particular, we compare the 
behaviours of multiqubit entanglement measures based (i) on the von Neumann entropy of marginal 
density matrices and (ii) on the linear entropy of those matrices. 
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I. INTRODUCTION 



The study of quantum entanglement is contributing both to the elucidation of the foun 
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dations of quantum mechanics and to the birth of new, revolutionary technologies 
A considerable amount of research has recently been devoted to the study of multiqubit 
entanglement measures defined as the sum of bipartite entanglement measures over all (or 
an appropriate family of) the possible bi-partitions of the full system 0, O, [7]. The aim of 
the present contribution is to explore some basic properties of highly entangled multiqubit 
states, and also of the "entanglement landscape" in their neighbourhoods. We compare the 
behaviours of two entanglement measures for multiqubit pure states, one based on the von 
Neumann entropy of marginal density matrices and the other based upon the linear entropy 
of those matrices. We also compare the performances of two searching algorithms for highly 
entangled states, based on different families of bi-partitions of the multiqubit system. 

II. MULTIQUBIT ENTANGLEMENT 

The genuine multipartite entanglement E of a A-qubit state can be expressed as 
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Here, E(i) stands for the entanglement associated with one, single bi-partition of the 
A-qubits system. The quantity E^ m ' gives the average entanglement between subsets of m 
qubits and the remaining N — m qubits constituting the system. The average is performed 
over the N^ art nonequivalent ways to do such bi-partitions, which are given by 



N Wart = ( ^ n ^ A/2, (3) 
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Different E^ represent different entanglement properties of the state. While E^ can 
attain its maximum value for a given state, E^ can be arbitrarily low for such state. This 



is why all these entanglement measures must be computed to capture all the entanglement 
properties of the state. The global multiqubit entanglement is given by the average of the 
[N/2] different £ (m) for any state |tf ). 

We will use two types on entanglement measures, El and E v n, respectively based on 
two different measures for the mixedness of the marginal density matrices pi associated with 
the bi-partitions: (i) the linear entropy Sl = 2 m™ i (1 — Tr[pj]), and (ii) the von Neumann 
entropy S v n = —Tr\pilogp/\. If one uses the linear entropy Sr. E L turns out to be the 
well known Meyer- Wallach multipartite entanglement measure 8] that Brennen |9j showed 
to coincide with the average of all the single-qubit linear entropies. This measure was later 
generalized by Scott [h]] to the case where all possible bi-partitions of the system where 
considered. 

The entanglement measure given in Eq. (1) is maximized by a state which has all its 
reduced density matrices maximally mixed. Although it is easy to verify that in the 3 qubit 
case \GHZ) complies with this requirement, the situation becomes much more complicated 
when systems with four or more qubits are considered. Higuchi and Sudbery proved that 
there is no 4 qubit state whose two qubit reduced density matrices are all maximally mixed. 
They conjectured that the 4 qubit state exhibiting the higher entanglement is 

\HS) = -j= |1100) + |0011) +w(|1001) + |0110)) +cj 2 (J1010) + |0101)) , (5) 
with ijj = —1 + Although it still remains unproven, several analytical 
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numerical [13| , Il4j evidences support the aforementioned conjecture. In the cases of 5 and 
6 qubits, states have been identified having all their reduced density matrices maximally 
mixed [la . . Finally, for 7 qubits there is numerical evidence suggesting that a recently 
discovered state is the one with maximal entanglement although, as in the 4 qubit case, no 



state of 7 qubits with all its reduced density matrices maximally mixed was found 



III. BEHAVIOUR OF E L AND E VN FOR HIGHLY ENTANGLED STATES OF 4 
QUBITS. 

The multipartite entanglement measures El based on the averaged linear entropies of 
the reduced density matrices are widely used, but sometimes it is more convenient to use 
an entanglement measure E v ^ based on the von Neumann entropy of the reduced density 
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FIG. 1: Probability Density Function for E v pj among those states that maximize El- 

matrices, although its computation is not as straightforward as the computation of El- Here 
we compare the behaviour of the entanglement measures El and E v ^ when searching highly 
entangled states of 4 qubits. Our results indicate that E vN is the best measure to use. As 
shown in Fig. 1, most states that maximize El are not maximally entangled according E v ^, 
even though they are all highly entangled states (most of these states have a value of E v n 
around 0.935). 

The study of the set of highly entangled 4 qubits states is of considerable interest because 
they represent the lowest dimensional system for which the non existence of the theoreti- 
cally maximally entangled state has been proved. Brown et al. [3] developed a numerical 
algorithm to search highly entangled states of multi-qubit systems and found a maximally 
entangled state of 5 qubits. However, when applied to 4 qubit systems their algorithm yielded 
a state (which we here call \BSSB4)) less entangled than the \HS) state previously discov- 
ered by Higuchi and Sudbery. This state is \BS3BA) = |(|0000) + |+011) + |1101) + |-110)), 
where |+> = ^(|0> + |1» and |-> = ^(|0> - 

A new and slightly different numerical algorithm was recently developed by us [14j that 
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FIG. 2: Entanglement values for those states in the neighbourhood of \HS) and \BSSB4) as a 
function of the overlap with them. 



has been successfully applied to find maximally entangled states in systems up to 7 qubits, 
including the 4 qubit \HS) state. Here we compare the behaviour of El and E v n as mul- 
tiqubit entanglement measures for highly entangled states of 4 qubits, through the study 
of the entanglement properties of the states living in the neighbourhoods of \BSSB4) and 
\HS). To such an end we first compute the average entanglement of states having given 
overlaps with \BSSBA). We considered, in total, a family of 15 000 000 states |^'), with 
0.95 < \(V'\BSSB4)\ < 1. A similar computation is done with a second family of 15 000 000 
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states |$') (this time close to \HS)) with 0.95 < \(&'\HS)\ < 1. The results are summa- 
rized in Fig. IIHI While both entanglement measures identify all the alluded states as highly 
entangled, El does not succeed in distinguishing the neighbours of \BSSB4) from the neigh- 
bours of \HS). On the other hand, the averaged E v ^ measure successfully distinguishes both 
families of states, and identifies the states in the neighbourhood of \HS) as more entangled 
than those related to the \BSSB4). Interestingly, the slopes of both curves depicted in the 
upper part of Fig. 2 (indicating the rate of decrease in entanglement as we consider states 
with decreasing overlaps with \BSSB4) or \HS)) are approximately the same. This suggests 
that the "entanglement landscapes" in the neighbourhoods of \BSSB4) or \HS) share some 
basic features. 

IV. ALTERNATIVE APPROACH TO THE NUMERICAL SEARCH ALGO- 
RITHM 

In Ref. [14]] we proposed a numerical search algorithm that was able to find maximally 
entangled states in systems up to 7 qubits, starting from an initial separable state. To 
find the maximally entangled state, the coefficients of the initial state are slightly modified 
to obtain a new one. The entanglement of the new state is computed, if it is larger than 
the entanglement of the previous state the new state is kept. Otherwise, the new state is 
rejected and a new, tentative state is generated. This iterative process is repeated until it 
converges to a maximally entangled state. At each iteration the entanglement given by Eq. 
(1) must be computed. Consequently, at each step we must evaluate as many E^ measures 
as non-equivalent bi-partitions the system has. This implies an exponential increase with the 
number of qubits of the computational resources needed to find the final state. Consequently, 
it is highly desirable to develop schemes to decrease the number of iterations needed to obtain 
the convergence and the time needed to perform each iteration. 



Number of qubits 


3 4 5 6 7 


E L 


1.0000 0.9445 1.0000 1.0000 0.9961 


E V N 


1.0000 0.9481 1.0000 1.0000 0.9948 



If a multiqubit state has highly mixed reduced density matrices corresponding to subsys- 
tems of [N/2] qubits, it is reasonable to expect that the same will happen with the reduced 
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density matrices describing smaller subsystems. For this reason, in order to optimize our al- 
gorithm, we have tried a modified scheme based on the maximization of E^J 2 ^ . The results 
of this experiment have been reasonably successful. For systems of 3, 4, 5, and 6 qubits 
the final highly entangled states obtained maximizing E v ^ are the same as those obtained 
maximizing E^j 2 \ This is a big improvement in our numerical algorithm, because in each 
iteration the number of bi-partitions to be considered is roughly reduced to the half, and the 
total number of iterations needed to reach the convergence are usually considerably reduced 
as well. For 7 qubits the E vN entanglement values of the states yielded by the i?j^ 2 '-based 
algorithm differ in the sixth decimal digit from the E entanglement value of the optimum 
state obtained maximizing E v ^. For 8 qubits the optimization algorithms, based either on 
balanced bi-partitions or on the global entanglement measure, do not converge always to 
the same state. The entanglement values (for different number of qubits) of the multiqubit 
states of highest entanglement considered in the present work are given in Table 1. 

V. CONCLUSIONS 

We have compared the behaviours of the multiqubit entanglement measures El and E v n 
based, respectively, on the linear and the von Neumann entropies. Our results indicate that 
E v n is better than El for the search of highly entangled states, because it discriminates 
between states that, while exhibiting the same value of El, have different degrees of entan- 
glement. We also found evidence that search algorithms based upon balanced bi-partitions 
are almost as efficient as those based on the complete set of bi-partitions. 
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